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The free energy of the Ginzburg-Landau theory satisfies a nonlinear functional differential equa- 
tion which is turned into a recursion relation. The latter is solved graphically order by order in the 
loop expansion to find all connected vacuum diagrams, and their corresponding weights. In this way 
^^ . we determine the connected vacuum diagrams and their weights up to four loops. 
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I. INTRODUCTION 



Recently, two of us (H.K. and B.V.d.B.) have determined the two-loop effective potential of 0(7V)-symmetric scalar 
quantum electrodynamics in 4 — e dimensions and its e-expansion B. Going to higher loop orders requires the 
^ ' calculation of many complicated Feynman integrals associated with vacuum diagrams. In this paper, we show how to 
find all these diagrams and their weights with the help of a simple recursive technique described in detail in Refs. {§-f|- 
The basics of the method were formulated in Refs. J7],||. For simplicity, the present work is restricted to the theory of 
a single complex scalar field coupled to electromagnetism, known as Ginzburg-Landau theory or as scalar quantum 
electrodynamics. In D = 3 dimensions, this describes the physics of superconductors. Without electromagnetism, we 
y—( ' recover the physics of superfluid helium. Here we shall restrict ourselves to the symmetric phase, where the O(N)- 
symmetry is unbroken, which describes the system above the critical temperature T c . The more general situation will 
be dealt with elsewhere. 

For a complex boson field, propagators are represented by oriented lines, as in ordinary quantum electrodynamics 
pj. When comparing the two expansions we must, however, drop a minus sign for each fermion loop. In addition, the 
Ginzburg-Landau theory contains four-point functions of the scalar field which are in principle treated in <fi -theory 
P|. However, since we deal here with complex fields, the weight of the corresponding graphs cannot be taken from 
Ref. p[. There is furthermore a completely new vertex in scalar QED: the quartic seagull coupling between two 
photon and two scalar fields. By a replacement of the photon lines by scalar lines, they become equivalent to the 
^^ ' diagrams coming from the <j> -vertex, apart, of course, from different weights. 

The paper is organized as follows. In Section || we define more precisely the theory to be studied. In Section [II 



wo 



introduce basic functional derivatives which allow in Section [V to derive a functional differential equation for the free 
energy. From this follows a recursion relation which allows to find all vacuum diagrams of a given order L from those 
of the previous orders. This equation is formulated graphically in Section M, where we also determine the vacuum 
diagrams and the corresponding weights up to four loops. 

II. GINZBURG-LANDAU THEORY 

The generalized Ginzburg-Landau theory deals with a self-interacting complex scalar field coupled minimally to an 
electromagnetic vector potential. The theory contains two coupling constants g and e. The field expectation value 
(</>) vanishes above T c and the theory has three types of vertices. 

The physics can be extracted from the partition function 

Z = jV(j)V^VAexp{-E), (1) 

with thermal fluctuations being governed by the energy functional 



E{<j,\ 0, A„) = t 1 G 12 1 2 + 1^123401020^4 + ^MD^A2 + eH- U2 Ai4<f> 2 + \ ■FwMM&fa- (2) 

Here and below, ovcrlined indices are used for the electromagnetic field. We work in a covariant gauge and assume 
the photon propagator to have the appropriate form (see Ref. [B). Using the same notation as in Refs. ||,§J, we keep 
all equations as compact as possible by assuming Einstein's summation convention not only for the internal degrees of 
freedom, but also for the space-time indices, for which repeated indices imply an overall integral sign. In this notation 
the functional matrices in the energy functional (0) have the following symmetries: 

D ^= D W> (3) 

Vl234 = V3214 = U1432 = V3412, (4) 

^T212 = -^2112 ■ (5) 



III. BASIC FUNCTIONAL DERIVATIVES 

The recursion relation will be derived by performing functional differentiations of the partition function with respect 
to the propagators and their inverse as well as the vertices. The basic properties of these derivatives are shown in the 
following equations: 
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Vertex derivatives 
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This is a direct extension of relations used in Refs. [|3|-|3|. The main difference comes from the fact that the previous 
scalar propagators were symmetric in the indices whereas the present propagators describing complex fields are not. 
By the chain rule of differentiations, we can then find the derivative of any functional with respect to the propagators 
and their inverse as well as the vertices. 

IV. FUNCTIONAL DIFFERENTIAL EQUATION FOR FREE ENERGY 

With the definitions given in the previous section, we are now prepared to derive the graphical recursion relation 
for the vacuum graphs. We start from the identity 



J 6<f>\ 



lexp(-E) =0, (17) 



and obtain 
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which, using the Ginzburg-Landau energy functional (0), leads to the linear equation 

jvcfrVtfVA U 12 - 4Gli<t> 3 ~ 1^3454^4^5 - \FmaAlAi4<l>a ~ ^lnMAn\ ex P (~ E ) = °- 
Rewritten in terms of functional derivatives with respect to the inverse of the propagators, this becomes 

(v<j>V<jtVA U 2 + G^-^t - f^ms^^T - e2F m3-^TT7^T + eH-^A,-^) cxp (-25) = 0. (20) 
J y dG 23 I dG 23 dG 45 6D— dG 23 dG 23 J 

The last part of the above equation may be replaced by a three-vertex derivative. However, proceeding in this way 
would not lead to an iterative generation of diagrams. For this, it is necessary to consider a second identity: 



IVcjyVtfVA— exp (-E) = 0, 



(21) 



from which we obtain 



jv<t>VtfVA [D^A- 2 + aJWitfa + ^FjSia A i4<f>a) exp (-E) = 0, (22) 

or, equivalently, 

jvfiVtfVA (dz}A 2 - e#ii2^JzT ~ eF ^JW~) exp ( ~ E) = °- 
Inserting the photon-field expectation value (p3|) into (pfj|), we obtain 
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This equation is linear, but it leads to a huge number of diagrams, both connected and disconnected. We remove the 
disconnected ones by introducing the free energy W as generating functional of the connected Green functions as 

Z = exp (W) , W = Wo + W int , (25) 

where Wq is the free field part 



Wo = -TrlnG" 1 - -TrlnD- 1 , 



(26) 



with Tr being a shorthand notation for the functional trace. Working with the series representation of the logarithm, 
we obtain directly the relations 
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where the index ordering is important for the complex scalar fields. 

Introducing the decomposition (Eq) in ( p4|) and using the relations (27)-(|29|), we obtain a nonlinear functional 
differential equation for the interacting part Wi n t: 
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With the help of Eqs. (0) and (p"3|), this equation becomes 
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(30) 
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From this functional differential equation, a graphical recursion relation can be derived. This is the subject of the 
next section. 



V. GRAPHICAL RECURSION RELATION 



From (pl|) , we can derive a graphical recursion relation for the connected vacuum diagrams. When considering the 
loop expansion of the interaction part W; n t , one term is of one loop number larger than the other terms, as we now 
show. 

The operators G12S/SG12 and Dj^S/SDj^ simply count the number of scalar lines I q and photon lines I 7 , respec- 
tively, in a given diagram. These numbers can be extracted from the number and type of vertices. Denoting by 
V g , V e and V e 2 the number of g, e and e vertices, we have the following counting rules. The Yukawa vertex -ffj 12 
has one photon line and two scalar lines, while the quartic photon-scalar vertex F^ 12 has two photon lines and two 
scalar lines. Furthermore two vertex lines are necessary to produce an internal line when combining vertices. We 
have the obvious relation 2/ 7 = 2V e 2 + V e . Taking the quartic scalar self-interaction V1234 into account, we have also 
2/ q = 2V e 2 + 2V e + 4V g . An odd number of photon fields gives no contribution to the free energy. Thus the vertex e 
enters with even power. The number of loops is then easily found to be L — 1 = V g + V e /2 + V e 2. Together, we have 
the counting rules 
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(32) 

(33) 

L-l = V* + ^+V g . (34) 

These relations can be inverted to give the number of each type of vertex as a function of the number of loops, scalar, 
and photon internal lines: 

V e 2 = 2(1 - 1) - J q , (35) 

^ = I q + I r -2(£-l), (36) 

Vg = L-l-/ 7 . (37) 

It is now clear that, to count the number of scalars in a given diagram, it is only necessary to know the loop order, 
as well as the number of the quartic e photon-scalar seagull vertices: 

I q = 2{L-l)-V e2 . (38) 

The vertices V e are not taken into account when counting the quartic vertices V e 2 , which, by definition, count the 
two photon- two scalar vertices only. The relation (pq) is interesting, since, for complicated diagrams, with a large 
loop-order, it is much less involved to count the number of quartic e photon-scalar seagull vertices than the number 
of scalar lines. 

We are now ready to demonstrate that equation (|3l]) allows for a recursive solution. We form the loop expansion 

oc 

Wtat^Xy-HfW (39) 

supposing that the vertices e and e are of order y/g and g, respectively. This is because the relevant parameter for 
the loop expansion is the inverse temperature (3 = l/(ksT) which appears as a coefficient in front of the energy. The 
inverse temperature is set equal to one in this work (see Eq. (||)). We could have restored it to show the loop counting. 
We would have seen readily that this would have been equivalent to the statement that e and e 2 are of order yfg and 
g, respectively. In Eq. (|39|), W^ L > is a sum over the different diagrams of a given order L: 

W(L) =^{-l) V3+V - 2 W {L4 \ (40) 

d 

where d distinguishes between different classes of diagrams of the same loop order. The factor (— l) v s+ v e 2 takes care 
of minus signs in the diagrams of the perturbation expansion. Since the number of Yukawa vertices is even, it does 
not enter this prefactor: (— l) v " — 1. Applying the scalar number operator on W; nt gives 

gJ-^ = f;^- i ^(-i)^^/ q (i,d)w^ ! (4i) 

which stresses that each class of diagrams, i.e., each topology, has its own scalar number. 

By performing the loop expansion (pS9|), the contributions W' > to the negative free energy consists of all connected 
vacuum diagrams constructed according to the following Feynman rules. A straight line and a wiggly line represent 
the free correlation of the scalar and the photon field, respectively: 

1^2 = G 12 , (42) 

T 2 = Djs. (43) 

The vertices are correspondingly pictured by 

x = - e L H ^^ (44) 
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X = -Z 2 L *T234, (45) 
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With these Feynman rules, the insertion of the loop expansion (B9|) into (pl|) leads to the following equation 
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which we have written directly as a graphical equation where the vertices already contain the coupling constants 
e,e ,g. The last four terms are the nonlinear ones. They start to give a contribution only at the three-loop level. 
(Although the overall summation starts at L' = 2, the nonlinear terms have an internal summation for L' € [2, L— 1], 
leading to a vanishing two-loop contribution.) The fact that only positive signs enter is a consequence of the Feynman 
rules ©-©. 

The term W ', corresponding to two-loop diagrams, is the initial condition to enter this equation. It needs not 
to be given because it is also contained in (|47|). Identifying the terms of order g, we see directly that only the first 
line survives. Integrating the equation gives the four diagrams in the right-hand-side, with the corresponding weights 
1/2, 1/2, 1/2, 1/2 obtained by dividing each graph by its number of scalar lines (compare Table I): 
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Equating powers of g, we end up with the following nonlinear graphical recursion relation for the vacuum diagrams: 
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+ E ^T=^ ,X 4 ^4 + E ^T^ a >-C 3 1— ( 49 ) 
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The recursive nature of this equation is due to the fact that the left-hand-side is of one higher loop order than the 
right-hand-side: To obtain the diagrams at order L + l, one applies the right-hand-side operators on the lower order 
diagrams. Since the left-hand-side contains the scalar number operator for each diagram, the weight of the diagram 
which have just been obtained has to be divided by the respective number of scalar lines, see Eq. (f4l|). We remind 
the reader that this number can be obtained knowing the number of seagull vertices of the respective diagram, see 

Eq - ©■ 

We have derived the vacuum diagrams of the Ginzburg-Landau theory up to four loops which involves the nonlinear 

part of Eq. (f49|). Indeed, the nonlinear terms only enter from the three- loop order. The resulting diagrams are given 

in Table 1. The effort needed to obtain them is considerably reduced compared to a calculation done with the help of 

external sources coupled linearily to the fields. Note that Eq. (MS) is suitable for an automatized symbolic computation 

which can be implemented as in Ref. H, such that one may proceed to higher orders without much effort except for 

computer time. 

A look at Table 1 shows that the diagrams and weights of the Yukawa part coincide with those from Ref. m , as it 
should. For the pure <fi -part, the diagrams are equivalent to those in Ref. Q, although the weights do not coincide, 
since we deal with complex scalar fields. 

Being in the possession of all Feynman diagrams we must still calculate the associated integrals in order to extract 
physical results. This will be done in a separate publication. In particular, we intend to compute the vacuum energy 
which determines the critical behavior of the heat capacity of a superconductor at the phase transition. 

VI. CONCLUSIONS 

In this paper, we set up a graphical recursion relation for obtaining the connected diagrams of the Ginzburg-Landau 
model which describes superconductors near the critical point. We have used our equation to obtain the diagrams up 
to the four-loop order. These diagrams will be needed to extend our two-loop calculations in Ref. H to higher orders. 
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TABLE I. Connected vacuum diagrams W > ni '™ 2 '"3J anc j their weights up to the four-loop order of the 0(2) 
Ginzburg-Landau model, wehre L denotes the loop order and m , ri2, n3 count the number of vertices V, -F, H, respectively. 
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